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A field-theoretical description of the behavior of a disordered Ising system with long-range inter- 
action is presented. The description is performed in the two-loop approximation in three dimensions 
using the Pade-Borel resummation technique. The renormalization group equations are analyzed, 
and the fixed points determining the critical behavior of the system are found. It is shown that the 
effect of frozen structural defects on a system with long-range interaction may cause a change in its 
critical behavior or smearing of the phase transition. 
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The effect of long-range interaction, which at long dis- 
tances is described by the power law r D<T , was stud- 
ied analytically in terms of the expansion [13] and nu- 
merically by the Monte Carlo method [46] for two- and 
one-dimensional systems. It was found that the effect 
of long-range interaction on the critical behavior of Ising 
systems is considerable for a < 2. The study [7] carried 
out for a three-dimensional space in the two-loop approx- 
imation confirmed the prediction of the e-expansion for 
a homogeneous system with long-range interaction. 

According to the results of [8, 9], the introduction of 
frozen impurities into a system leads to a change in its 
critical behavior. In this connection, it is of interest to 
determine the effect of frozen structural defects on the 
critical behavior of a system with long-range interaction. 

This paper describes the critical behavior of a three- 
dimensional disordered system a long-range interaction 
for different values of the parameter a. 

The Hamiltonian of a system with long-range interac- 
tion has the form 

Ha = \J dDq{T ° + q ^ S « S -i + \j d D q^T q S q S- q 

(1) 

where S q denotes the fluctuations of the order parameter, 
D is the space dimension, tq ~ |T — T c \, T c is the critical 
temperature, uq is a positive constant, and Ar q is the 
random field of impurities of the random temperature 
type. 

The critical behavior strongly depends on the param- 
eter (7, which determines the rate of interaction decrease 
with distance. According to [1], the effect of long-range 
interaction on the system behavior is considerable for 
< a < 2, while, for a > 2, the critical behavior is 
equivalent to that of a system with short-range interac- 
tion. Therefore, only the case of < a < 2 is considered 
below. 



For a low impurity concentration, the distribution of 
le random field Ar q can 1 
described by the function 



the random field Ar q can be considered as Gaussian and 



(2) 



where A is the normalization factor and 5q is a posi- 
tive constant proportional to the concentration of frozen 
structural defects. 

Applying the replica procedure for averaging over ran- 
dom fields caused by the frozen structural defects, we 
obtain the effective Hamiltonian of the system 
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Hr = - d D q{T + q°)Y,StS- q (3) 
J b=i 

S m f 
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The properties of the initial system can be obtained in 
the limit of the number of replicas (transforms) m — * 0. 

Applying the standard renormalization group proce- 
dure based on the Feynman diagram technique [10] with 
the propagator G(k) = 1/(t + \k\ a ), we arrive at the 
expressions for the functions (3 u ,f3si lip and 7t, which de- 
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termine the differential renormalization group equation: 



0u 



= -(2a -D)u 1 - 36u J + 245 J 

+ 1728 (2 Ji - J 2 - ^Gju 2 - 
- 2304(2 Ji - J 2 - \g)u5 + 672(2 Ji - J 2 - ^G)c5 5 



ft = -{2a-D)S 



1 - 24uJ + 168 J 



+576(2 J! - J 2 - -Gi)w 2 - 
- 1152(2 Jx - J 2 - ^G>5 + 352(2Jx - J 2 



the cutoff parameter A and considering the ratios Ji/Jq 
and G/Jq we obtain finite expressions in the limit A — > 
(4) 00. 

The values of the integrals were determined numeri- 
cally. For the case a < D/2, a sequence of the values 
of Jx/Jq and G/Jq was constructed for different A and 
approximated to infinity. 

The critical behavior is completely determined by the 
stable fixed points of the renormalization group trans- 
formation. These points can be found by equating the 
/3-functions to zero: 
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/JiK>«2) = 0, fo(v* 1 ,v* 2 )=0. 



(7) 



7t = (2a -D) - 12?iJ + 4<5J + 288^2 Ji - J 2 - ^G^u* 

- 288(2 Jx - J 2 - ^G)m(5 + 32(2J X - J 2 - ^G)S 2 
7 V = (2cr- J D)64G(3u 2 -3it(5 + (5 2 ), 



The stability requirement for a fixed point reduces to the 
condition that the eigenvalues bi of the matrix 



dvj 



(8) 
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(i + kr) 2 a + br) 
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'(i + i g 2 +p 2 + 2p,7r /2 )' 

f d D q 

J (1 + kr) 2 ' 

9 f d D qd D p 



d\k\° J (1 + \q 2 + k 2 + 2kq\°) 



(l + \p\ a )(l + \q 2 +p 2 + 2pq}°/ 2 ) 

Let us redetermine the effective interaction vertex: 

vi = u ■ J , v 2 = 5 • Jo- (5) 

As a result, we obtain the following expressions for the 
functions (3, -f v and 74: 



A = -(2a - D) 



+ 1728(^2Jx - 1 



1 - 36^x + 2Av 2 

1 % 



(6) 



2304(2Jx - 1 - -G)v lV2 + 672(2 Ji - 1 - -G)v 2 



fo = -(2a-D)8 



1 - 2Avi + 16u 2 



+576(2Jx - 1 - -G)vl 

-1152(2Ji - 1 - \g) Vi v 2 + 352(2J7 - 1 



lie in the right half-plane of the complex plane. 

The critical exponent v characterizing the growth of 
the correlation radius in the vicinity of the critical point 
((R c <~ \T — T c \~ v )) is obtained from the relation 

u = -(l + lt )'\ 
a 

The Fisher exponent 77 describing the behavior of the 
correlation function in the vicinity of the critical point in 
the wave-vector space ((G ~ k 2+v )) is determined from 
the scaling function 7^: rj = 2—a+-/ v (vl, v 2 ). The values 
of other critical exponents can be determined from the 
scaling relations. 

It is well known that the perturbation series expan- 
sions are asymptotic, and the interaction vertices of the 
order parameter fluctuations in the fluctuation region are 
sufficiently large, so that Eqs. (7) applies. Therefore, to 
extract the desired physical information from the expres- 
sions derived above, the PadBorcl method generalized to 
the many-parameter case was used. The corresponding 
direct and inverse Borel transformations have the form 

00 

f(vi,v 2 ) = £ c iui2 v^v^ = J e- t F(v 1 t,v 2 t)dt, 

*1,*2 (9) 

For the analytic continuation of the Borel transform of 
a function, a series expansion in powers of the auxiliary 
variable is introduced: 



7t = (2a-D)\- 12 V1 + 4v 2 + 288 f 2 J x - 1 - -G)v 2 F( Vl ,v 2 ,6) = ^0* ^ -^^ 1 ^M il+i2;fe , (10) 

L V 3 / 1 — n , .■ K - 



-288(2Ji - 1 - \g)v iV2 + 32(2Ji - 1 - )-G)vl 

7 y = (2a - D)64G(3vl - 3 Vl v 2 + v 2 2 ). 

Such a redetermination makes sense for a < D/2. In this 
case, Jq, Ji and G are divergent functions. Introducing 



k— £1,12 

and the [L/M] Pade approximation is applied to this se- 
ries at the point 9 = 1. This approach was proposed and 
tested in [11] for describing the critical behavior of sys- 
tems characterized by several interaction vertices of the 
order parameter fluctuations. The property (revealed in 
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a 


«i 


V2 


bi 


fc 2 


1.5 


-0.395432 


0.951135 


84.530 


59.517 


1.6 


-0.227628 


0.594810 


45.302 


32.575 


1.7 


-0.045234 


0.274890 


13.235 


3.915 


1.8 


0.064189 


0.046878 


0.626* 


0.626* 


1.9 


0.066557 


0.040818 


0.559* 


0.559* 



[11-14]) that the system retains its symmetry when us- 
ing the Pade approximants in the variable 6 is essential 
in the description of multivertex models. 

The table shows the stable fixed points of the renor- 
malization group transformation and the eigenvalues of 
the stability matrix at a fixed point for the parameter 
values 1.5 < a < 1.9 One can see that stable fixed points 
exist in the physical region (vl,v% > 0) only when the 
long-range interaction parameter is a a > 1.8. The cal- 



culations showed that, for any a < 1.8, the stable points 
of a three-dimensional impurity system are characterized 
by a negative value of the vertex. 

The calculation of the critical exponents provided the 
following values: 

o-=1.9, v = 0.706786, r/ = 0.1344048 (11) 
ct = 1.8, v = 0.732789, r\ = 0.250978 

Thus, for a three-dimensional Ising system with long- 
range interaction, the introduction of frozen structural 
defects leads to a change in its critical behavior if the 
long-range interaction parameter is a > 1.8 and to smear- 
ing the phase transition if the parameter is a < 1.8. 

The work is supported by Russian Foundation for Ba- 
sic Research N 04-02-16002. 
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